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Abstract 

We calculate the excursion and meander area distributions of the elastic Brownian 
motion by using the self adjoint extension of the Hamiltonian of the free quantum particle 
on the half line. We also give some comments on the area of the Brownian motion 
bridge on the real line with the origin removed. We will stress on the power of self 
adjoint extension to investigate different possible boundary conditions for the stochastic 
processes. We discuss also some possible physical applications. 
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1 Introduction 

In this paper we study some area distributions of elastic Brownian motion [Ij , such as elas- 
tic Brownian excursion and elastic Brownian meander. The goal apart from the calculation 
of the area distribution of elastic Brownian motion is to give a unified framework to study 
different area distributions for Brownian motion with generic boundary conditions. This 
generalization is in close connection with the concept of the self adjoint extension in quan- 
tum mechanics. The self adjoint extension gives a reasonable classification for the possible 
boundary conditions of the quantum particle and so for the stochastic process. 
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Different area distributions of Brownian motion in one dimension were calculated by the 
mathematicians in the last century. The Brownian motion area was first calculated by Kac 
[2]. The Brownian excursion area was calculated by Darling [3j and Louchard [5 and the 
Brownian meander area was calculated by Perman and Wellner [5] . For the applications and 
other distributions of the areas see [H [7] and also the complete review by Janson [8] , and 
references therein. 

The quantum mechanics methods first used by Comtet and Majumdar [9j to re-drive the 
different area distributions of Brownian motion. This method which is more common in the 
physics literature apart from simplicity can give a unique way to calculate the different area 
distributions. It is also very useful to generalize the results to more complicated stochastic 
processes. The self adjoint extension of the Hamiltonian operator gives us the variety of 
different possible boundary conditions in the presence of the boundary [10]. In this letter 
we will focus on the area distribution of the elastic Brownian motion which is equivalent to 
the free quantum particle on the half line. We also calculate the area distribution of the 
stochastic process equivalent to the free particle on the real line with a hole at the origin. 
Of course the area distribution is just one of the myriad of possible distributions that one 
can investigate for the elastic Brownian motion and the Brownian motion on the pointed real 
line. We will summarize some of these distributions in the last section. 

The paper is organized as follows: in the second section we define the elastic Brownian 
motion and its connection to the quantum particle on the half line. In the third section we use 
the same method as [9j to calculate the area distribution of the elastic Brownian excursion 
and meander. Different limits of our calculation will give the well known results. We will 
conclude this section with some immediate application of our results and giving some hints 
about other useful Brownian functionals. In the forth section using the results of the section 
three we will give the distribution of the area for Brownian bridge with a point defect at the 
origin. Finally in the last section we will summarize our results and possible future directions. 

2 Elastic Brownian motion and quantum mechanics 

The elastic Brownian motion is the natural generalization of the Brownian motion in the 
half line with special boundary condition on the origin. To define the process we need first 
to introduce the local time. The definition of the local time of the path oo at the point a is 
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as follows 



1 , 1 



T 



ti{a) := -lim- / l,+,(5,)o!s, (2.1) 

where lx+e{Bs) is the indicator for the time that the process is in the interval [0,x + e]. 

One can naively write the above equation as an integral over a delta function as ti{uj,0) = 

6{x{t))dt. The local time has dimensions of the inverse of velocity and can be written as 

ti{a) = for the discrete random walk, where n is the number of times the path hits the 

origin. Then one can write the exponential of the local time as 

—Itt Itt 
exp( Axf. (2.2) 

J] T] 

To go to the discrete level we multiplied ti with that comes from the central limit 

theorem. The equation (12. 2p means that by considering 1 — "^ti as the probability for a 
single reflection from the origin it is possible to interpret the exponential of local time as 
the probability that a particle on a given path is reflected from the origirQj. Then the Green 
function of the elastic Brownian motion is just by the following expectation value 

-27r 

G(x,y,r) =<exp( ti) > . (2.3) 

The elastic Brownian motion is in close connection with the quantum particle on the half line. 
The Hamiltonian operator of the quantum particle on the half line has self adjoint extension 
with the following boundary condition 

m - (2.4, 

The above boundary condition is called Robin boundary condition. The energy of the particle 
is E]. = and the wave functions are 

V'fe = '^cosikx + 6k), (2.5) 

where tan((5fe) = ^ and 5^ is the phase shift corresponding to the s-wave. The solutions are 
normalizable and complete. 

The Hamiltonian is self adjoint for all of the real values of rj but to avoid the cases with 
bound states we will consider just non-positive ij. t/ = is the Dirichlet boundary condition 



^The name elastic Brownian motion comes from this property of tiie process. 77 = oo is the reflecting 
barrier, the particle will be reflected from the boundary with probability one. = is the absorbing barrier, 
the particle will be absorbed after hitting the boundary. Other cases between two extreme cases called elastic 
barrier. 
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and — > — cxD is the Neumann boundary condition. The Green function with respect to the 
solutions of the Hamiltonian has the fohowing form 

/•oo 

G^{x,y,t)= dke-"''''ip{y)^*{x). (2.6) 
Jo 

Using the above equation the Green function has the following form for arbitrary values of 
the self adjoint extension [TTl [12l [13] 

Gr^{x,y,t) = GF{x-y,t) + GF{x + y,t) + — [ dwe^'"GF{x + y + w) r? < 0; (2.7) 

^ Jo 

Grj{x,y,t) = Gf{x - y,t) + Gf{x + y,t) / dwe ^ Gf{x + y - w) 

V Jo 

+ ^![e'^e--("+^) rj>0; (2.8) 
V 

GF{x-y,t) = ^=e*(^-s')'/2t. (2.9) 
For the special cases, Dirichlet and Neumann the results are as follows 

Grj=o{x, y, t) = Gf{x -y,t)- Gf{x + y, t); (2.10) 
G^^_oo(x, y, t) = Gf{x -y,t) + Gf{x + y, t). (2.11) 

One can use the above equations to get the Green function of the elastic Brownian motion by 
just Wick rotation. The important point of this section is the possibility of using quantum 
mechanics language to describe the elastic Brownian motion. The other interesting point is 
the possibility of extending this equality in to the level of path integral representation [111 [T2| 
I13j . In the next section we use this correspondence to calculate different area distributions 
of the elastic Brownian motion. 

3 Area of the elastic Brownian motion 



In this section we will solve the problem of the area for the restricted Brownian motion, 
in particular we will solve the problem of the area distribution of elastic Brownian excursion 
and elastic Brownian meander. In the extreme limits we will get the well known results. 

3.1 The Area Under the elastic Brownian excursion 

In this subsection we calculate excursion area of the elastic Brownian motion. The definition 
of the excursion area is as follows: take an elastic Brownian process starts at x(0) = e 
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and comes back after time T to x{T) = e, without crossing the origin between. We are 
interested to the probabihty density P(yl,T, e) of the area A = x{T)dT for a fixed e and 
then finally take the limit e ^ 0, it plays the role of the regulator and can be treated 
independent of the limit in the local time process. This regularization is just necessary for 
the Dirichlet boundary condition but we are happy to keep it in our calculation to see its 
relevance in the calculation of Dirichlet boundary condition. To do this calculation we map 
the problem to the quantum mechanical problem. Now we use the method of Comtet and 
Majumdar [9j to calculate the excursion area distribution. To satisfy the constraint that 
process stays non-negative between and T one can multiply the above measure with the 
indicator function nr=o ^[■^(^)] ^-^P (~f^^') '^^lich 6 is the step function. The distribution 
P{A, T) of the area under the elastic Brownian excursion can then be expressed as the 
following quantum mechanical problem 

T) = ^ / Vx{T)e- ^« TT e[x{T)]5{ / x(r)dr - A), (3.1) 

Jx{0)=e ^J-Q Jo 

where is the normalization and corresponds to the quantum mechanics of a particle with 
infinite wall at the origin and zero potential at the positive real line 

Zl=<e\ e-^°^ I e > . (3.2) 

The Hamiltonian, Hq is equal to the self adjoint Hamiltonian that we discussed in the previous 
section. After integration, for small e we have 

^ 2(^ - 2vr6)2(— - 4e'^Erfc(^^)). (3.3) 

The integral for the Dirichlet and Neumann cases are 

Z% - -^e'T-'/' + Oie')- = 0, 

V 27r 



Z^ ^ ^:^ + 0{e'y, 7?^-oo. (3.4) 

The above partition functions are just the probability that an elastic Brownian motion goes 
from x{0) = e to x(r) = e in time T without crossing the origin. Taking the Laplace transform 
P(A, T) = P{A, T)e-^'^dA of both sides of the the equation gives 

Ze Jx(0)=e f-J-Q 
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In the numerator we have the propagator < e | e~^^'^ \ e > where Hi = — ^(^)^ + V{x) 
and V{x) = Xx for x > and infinite for x < 0. We absorb the Dirac delta function in to 
the boundary condition as the case of quantum particle on the half line, in the other words 
we consider the self adjoint extension of this operator. The boundary condition of the wave 
function after self adjoint extension is the same as (j2.4p . The solution of the Schrodinger 
equation is given by the Airy function as follows 

.*.'(.) = . ■^'((^^)'^^(^-^/^)) (3.6) 

^S„°"AiHi2xy'Hy-E/x))d, 

Using the boundary conditions one can determine the discrete eigenvalues as follows 

_ 2-1/3 .2/3 . M-4) _ .3 7^ 

Unfortunately since the above equation is transcendental the exact form of for the generic 
boundary condition is not available. However, for the Dirichlet and Neuman boundary con- 
ditions the solutions are just the magnitude of the zeros of Ai{z) and Ai'{z) on the negative 
real axes respectively, we show them by — and — c?°. The first few real roots of Ai{z) are 
approximately -2.33811, -4.08795, -5.52056 and -6.78671, etc. The first few real roots of 
Ai'{z) are approximately -1.018792, - 3.248197, -4.820099 and -6.163307, etc. Then the 
wave functions in these two cases are 

. ,3.S) 
^ (2A)V. ^^^-^) . (3.) 

To get the above results we used the identity Ai^{x)dx = —xAi'^[x) + Ai''^{x). Using the 
energy eigenvalues and wave functions one can write the equation (|3.5p as follows 



oo 



P^iX,T) = = 1^(^)1'^-" ' (3-10) 
For the Dirichlet and Neumann cases after considering small e we have 

oo 

P\\T) = ^Ar3/2^e-2-'''^'''^?^, (3.11) 

i=l 

P-iX,T) = #(2A)V3f;i_e--^^^^^^^'=rT. (3.12) 
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One can write the equations (j3.1ip and (3.12) in the scahng form as follows 

oo 

pO(A,r) = 5^/2^^6-2"^/'^'^''=", (3.13) 

i=l 

P-(A,r) = ny|2V3f;_Le-^-^/^«^-, (3.14) 

i=l * 

where s = XT^/'^ and u = T^/^A^/'^. It is possible to do the inverse Laplace transform of the 
functions (j3.1ip and (3.12) explicitly and find the distribution of the excursion area. To do 
so we need the formula 

£"i[exp(-.A");yl] = -^M„(.A-'^) (3.15) 
where {z) is the well known Wright function given by the following series 

2 

M ( \ 1 (-^)'^,, (fc + l) , . ,, (fc + l) , ^ 

fc=0 

For the Dirichlet boundary condition the inverse Laplace transform for T = 1 gives 

P\A) = S,dA{A-\M{%A-l,i^) = 

^ t{ 23 3 

( j ty( j. (3.17) 

^ fc=i 

where U{a,b,z) is the confluent hypergeometric function. 

To calculate the moments of the area one can work in the Laplace space and then use the 
equality r(n) < A"" >= P(A, T)A"-MA. To do the calculat ions we need to first define 
the generalized Riemann function A'^{s) = X^i^i t4t7 where comes from the equation (j3.7p . 
The above relations help us to calculate the moments explicitly as follows 

< A" V2^2^ "^''"^> A°(Hz!i±I>). ,3.8, 

1 (2 — n) 2 

It was shown in [14j that the moments after regularization are 

4-n r(n + 1 

r(^) 

where Kn is by the following recursion relation 



<A->=V2^2'^'^^^K^ (3.19) 



3 — 4 

Kn = -^^Kn-l + KjKn-j, s > 1, (3.20) 
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the first few values are Kq = — ^, Ki = |, K2 = ^ and K3 = j^. Then the first few 
moments of the excursion area are 



,n .1 \/27r ,9 5 .o 15v2vr 

<A°>=1, <^i>=X_, <A^>=-, <A^>=^^,.... (3.21) 

There is also a nice relation between the Airy zeta function and Kn as follows 

A«(^) = -l^^A-(„). (3.22) 

^ 2 ^ 3 sin(7rn) ^ ^ ^ ^ 

For example we have the following limits 

limnAO(|(n-l)) = A, A^O) = ^ (3.23) 
For the Neumann boundary condition we need another Laplace transform pair 

£-^[X-''exp{-sX-'');A] = ^^(f>{a,a;-sA''); -1 < a < 0, s > 0, < a < 1, (3.24) 

where i?!)(a, a; z) = X^fcLo k\v{ak+oi) generalized Wright function defined for the a > — 1. 

Using the above formula the area distribution has the following form 

^°°(^) = 2V3f;_L^(_L_^(^,|_2-V3,c«^^^)). (3.25) 

i=l * 



The moments of the area for T = 1 can be written as 

< A" >= ":;r;. ' > '-(:— :^). (3.26) 



3^r(l^) 3-3n, 
2^ r(-n) ^ 2 



These moments are the same as the moments of the Brownian bridge and can be regularized 
in the same way [8j . Then the moments are 

^2-"/2r(l + n) 
2 ^ 

where Dn is by the following recursion relation 



< ^" >= r(^^ — - °' ^^-^^^ 



3?7- — 2 2^ " 1 

i=i 

The first few values are Do = 1) Di = ^, D2 = ^ and = Then the first few moments 
of the area are 
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It is not possible to find exact probability distribution for the generic r] because for Robin 

boundary condition c'- is not independent of A and since they are related by non-algebraic 

relation it is not possible to get with respect to A explicitly. However one can follow the 

calculations in the level of Laplace space. The wave function has the following form 

, (2A)VMi((2A)V3a; - c;') , , 

^clAi^{-cl) + Ai'\-cl) 

The Laplace transform of the distribution of the area is 

To pursuit the calculation let us consider small r/s. One can write = + dcj where 
is the small perturbation around the zeros of the Airy function. After expansion of (j3.7p the 
perturbation is 

(3.32) 

For small r/ one can also expand Ze as follows 

Zl^:^,+0{v'). (3.33) 
Then the P(A, T) after expansion is 

oo 

P\\T) ^ V2^XtIJ2^~'"'''''''''^~'^'- (3-34) 

i=l 

The inverse Laplace transform of the function after using the equation (|3.15p is 

92/3 °° / r, 5 n 2 2 \ 

P^{A) ^ V2^— ^ c^dA {A - ±)-iM{MA - |-)-^ 3) • (3.35) 
k=i ^ 23 / 

The moments of the area can be find by the same method as before by just replacing A in 

the equation (IXT9|) hy A - ^. 

For the large rjs the same calculation can be done. The partition function of the elastic 

Brownian motion for large rj is 

Zl^\f^-- + Oi^). (3.36) 
y -iri r] r/^ 



The perturbation of c- after expansion of the equation ()3.7p is 

—Itt 

6dl ^ ^^TTT- (3-37) 

' c°°r/(2A)i/3 ^ ^ 

Unfortunately since the perturbation of the energy is dependent on the energy level we are 
not able to find simple equation for the distribution of the area in this case. 



3.2 The area under the elastic Brownian meander 

The definition of the elastic Brownian meander is similar to the elastic Brownian excursion, 
the only difference is for the elastic Brownian meander we do not need to force the process 
to come back to the starting point. In this case the partition function is 

/•oo 

= I db< 6|e-^o^|e > . (3.38) 

One can show that 







,1 —R„Ti I ^ I J- 1 (g+fc)^ \ 

<6|e^''-'|e> = W^-^ ( e"2^z^ + e"2^z^ j + 

f " Erfc[{^)'/^-^ Z -1/2. 3.39 

For small rj one can expand the error function up to the second order 

Taking the first order correction with respect to the and integrating over b gives 

Zl^Erfi^) ^e-^. (3.41) 

For r/ ^ oo it is easy to get 

^ 1- (3.42) 

Similar to the calculation that we did in the elastic Brownian excursion case one can write 
the Laplace transform of the distribution of 

1 

pV{^X,T) = -rr db <b\e-^'^\e> . (3.43) 

Jo 

Using the wave function (I3.30p one can get 

After expansion of the function with respect to the e and r] the first correction appears in the 
spectrum as follows 



P^(A, T) = 0F2-i/6(Ar3/2)i/3 ^ Bic^)e~^-'''^''''^'^. (3.45) 

1=1 
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^ ^0 -^■'^ 71 C\ Y) 

where B[c-j = ^^,^_^o) — and c- = c- + 6c- . The distribution of the area after inverse 
Laplace transform is 

P^{A, T) = V^2-V6tV2 g B{c^)^{ \ <t>{^, I -2-V3,.r(^ _ |-)^))(3.46) 

1=1 ^ 27r'' 

The continuation of calculation is now straightforward we just need to use the well known 
results for the Brownian meander. The moments of the area for Brownian meander, i.e. 
7/ = 0, for r = 1 is 

< A- >= V^2'-/'^^^Qn. (3.47) 
Qn satisfies the following recursion relations 

n 

i=i 

/3n = a„ + ^(2n-l)/?„_i, (3.48) 
6-2" r(3n + l/2) 

Oin — 



r(n + 1) r(n + 1/2) ■ 

The first few values are Qo = 1; Qi = Q2 = f| and = Then the first few values of 
the moments of the area are 

,n .1 2) pK ,9 59 ,n 465 pK , . . 

To get the results for the small rj one needs to replace A with j4 — ^ in the (j3.49p . 

The next interesting case is the Neumann boundary with r/ — > 00. The Laplace transform 

of the distribution of the area after a little algebra is 

00 

i=l 

where ki = ePQAtXc"") ^^^^ ^^{^) — ^i{y)(iy- This distribution is exactly the same as the 
distribution of the area of the Brownian motion, i.e. Jq \Bt\dt. This is not surprising because 
the absolute value of Brownian motion is like considering the area in the presence of totally 
reflecting boundary condition. Using the inverse Laplace transform one can get 

P^i^^T) = ^4H7r E-^-rM|(— ^) (3.51) 

1=1 

Using the well known results for the area of the absolute value of the Brownian motion [8] 
the moments of the area can be written as 

2-"/2r(l + n^ 

T(l^ 
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where L„ satisfies the following recursion relation 



n 



6i + i 

6j-l 




(3.53) 



i=l 



and CKj and I3j are the same as in the equation (j3.48p . The first few values are Lq 
-L2 = I and L3 = ^ . Then the first few values of the moments of the area are 



1, Li = 1 




(3.54) 



3.3 Some applications 

Different applications of the area distributions of Brownian motion in computer science, 
graph theory, fluctuating one-dimensional interfaces and localization in electronic systems 
were already discussed in length in many papers, see |16^ [T7] and references therein. The 
Airy distribution function and its derivative appear extensively in all applications. In this 
subsection we will summarize some immediate simple applications of our extended Airy dis- 
tribution. We will discuss the distribution of the average distance of a particle from a disorder 
with point interaction in three dimensions and the distribution of the average relative height 
distribution of interacting interfaces in two dimensions. We will also discuss one more elastic 
Brownian functional distribution related to vicious walkers interacting with the boundary. 

The first immediate application comes from the equality of the norm of the 3-dimensional 
Brownian motion (called three dimensional Bessel process) and Brownian motion on the half 
line [U |10] . This is very easy to show by considering the radial part of the Fokker-Planck 
equation of three dimensional Brownian motion. Now consider a Brownian motion moving 
in three dimensions in the presence of the disorder at the origin , one can map the system 
to the problem of one particle on the half line. The most generic point interaction between 
disorder and particle comes from the self adjoint extension of the Hamiltonian of the free 
particle in punctured three dimensional space [10] which is equal to the free particle on the 
half line with the boundary condition that we discussed in section 2. Then it is easy to see 
that the area distribution that we calculated is just the average distance distribution in the 
period T between the disorder and the free particle with the generic point interaction. 

Another simple application, which is in the close connection with the previous example, 
is the interacting interfaces. A path of Brownian motion in the x — t space is just like an 
interface with zero roughness exponent. One can also look at this interface as an ensemble of 
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growth models such as Edwards- Wilkinson model. Consider now two non-crossing interfaces 
in the region [0, L] with the similar boundary conditions. This problem is equivalent to the 
problem of two free particles on the real line. Consider, like the previous example, point 
interactions between the particles. The interaction between the particles is equivalent to the 
interaction between the interfaces. Then the area distribution that we calculated for the 
elastic Brownian excursion is just the distribution of the average distance of the interfaces in 
the interval [0, L\. One can also relax the boundary condition in one of end points and consider 
different boundary conditions for the different interfaces and use the results corresponding to 
elastic Brownian meander area. 

Since the elastic Brownian motion is the generalization of Brownian motion in the presence 
of the boundary we believe that all the possible applications should deal with the boundary 
interaction or point interaction between two particles. In this paper we just calculated one 
possible functional of the elastic Brownian motion, the area. However, there are many other 
functionals that can have interesting physical applications in the study of the interacting 
non-crossing walkers or interacting interfaces. We will discuss some of these functionals in 
the conclusion of the paper and give here one more example with more detail. 

Consider the problem of p non-crossing walkers, which has application in describing do- 
main walls of elementary topological excitations in the commensurate adsorbed phases close 
to the commensurate-incommensurate transition [18j, in the presence of the boundary. One 
interesting quantity is the maximal height distribution of the top walker that was already 
calculated exactly for the vicious walkers in [19] . Walkers are vicious if they annihilate each 
other when they meet. For simplicity we will discuss the simplest case p = 1. Consider H 
as the maximal height of the walker in [0, 1] then one can define the cumulative distribution 
as P{M) = Prob[H < M]. We define N{e,M) as the probability that the walker do an 
excursion in the period [0, 1] starting at e and coming back to the same point staying within 
the interval [0,M]. The cut-off e is just necessary for the Dirichlet boundary condition as we 
discussed before. Since for this case the result is already known [20] we will put e = here- 
after. Then it is easy to show that N{0,M) = J2e \i^E?e.-^. The wave function is as (12.51 ) , 
i.e. 'i/'E = \j^cos{kx + 6^), with = and k is the solution of the equation cot(A;M) = ||. 
Finally one can write the cumulative distribution as 



where comes from the equation (|3.3p . For small rj and small M with M < rj one can 
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simplify the equation as P{M) « ^^kQe~''o/'^ ^ wliere ko — — — . 

Generalization of the above results to arbitrary p is straightforward, however, it is rather 
cumbersome. From the perspective of our study in this paper one can generalize this problem 
in two directions: the first possibility is considering non- vicious walkers, interacting domain- 
walls, with a Dirichlet boundary condition on the wall. The second possibility is considering 
vicious walkers with non-trivial interaction with the boundary. 

4 The area of the Brownian bridge on the Une with a point 
defect 



(4.1) 



In this section we will summarize some results come from the area calculation for the 
Brownian bridge on the line with a point defect. The quantum mechanical counterpart was 
discussed lengthly in the literature and it is called the free particle on a line with a hole [15j. 
The functional integral for this problem was discussed in [13] and it is based on the different 
local times of the particle in the two sides of the origin. The most general boundary condition 
that respects the time reversal symmetry for the defect on the origin, after using self adjoint 
extension theory, is 

^+(0) / \-S -7/ V V'-(O) 
with a constraint 07 — (36 = 1. For simplicity we will consider some special cases. It is easy 
to see that for 5 ^ cxd the boundary condition decouples and we will have two decoupled half 
lines and so the area distribution is as the previous section. 

The next simple case is a = 7 = — 1 and 6 = which is equal to the free particle on the 
line with the following delta function potential 

V(x) = -^6{x). (4.2) 

To avoid the bound state we consider non-positive /?. To calculate the area of this kind of 
Brownian bridge one can use the method of previous section. Interestingly the results are 
very similar to the previous section. The energy of the particle is = ^k"^ and the wave 
functions are 

'^k{x) = ^J^cos{k\x\ + 6k), (4.3) 
14 



where tan{5k) = Then is the same as (|3.3p after replacing rj with The wave 
functions in the presence of the potential A|x| can have different parities, they are 

i^iix) = / ,-0, ^ , 4 .,9/ even parity, (4.4) 

y'C;3jA«^(-C;3j) + Al'''{-Cj3i) 

M^) = sgn(x)2-V2(2A)V6M^^)^^!N^fz2£i), odd parity, (4.5) 

where is the same as after replacing r/ with Odd wave functions do not contribute 
in the distribution of the area. The result for the even part is the same as the result for the 
elastic Brownian excursion. It is easy to see that /? = is like the free Brownian motion and 
so the distribution of the area is like Brownian bridge or like elastic Brownian excursion with 
Neumann boundary condition, i.e. (j3.25p . The P = — oo is like the the Dirichlet boundary 
condition and the distribution of the area is (|3.17p after considering rj = 0. 

The other simple boundary condition comes from a = 7 = — 1,/3 = and 5 7^ 0, in the 
other words 

V';(0)-V'-(0) = 0, (4.6) 
V'+(0)-V-(0) = -6i^'_{0)- (4.7) 

The energy of the particle is = ^k"^ and the wave functions are 

ipkix) = sgn{x)^J^cos{k\x\ + 6k), (4.8) 

where tan((5fc) = Then is the same as (|3.3p after replacing rj with 7r6. 

The wave functions in the presence of the potential A|x| can have different parities, they 

are 

,5/ ^ , ,2-V2(2A)l/6^^((2A)V3|^| 

ipf{x) = sgn(x) ^ ^ u y I I parity, (4.9) 

JcfAi^{-cf) + Ai'^{-cf) 

Mx) = (2A)^/^ ^^^2tcT^ ' even parity. (4.10) 

The above calculation shows that in this case the distribution of the area can be calculated 
by adding two terms. In this case the even parity has contribution for the distribution of 
the area which is equal to the Neumann boundary condition. Replacing rj with irS in the 
formula of elastic Brownian excursion will give the contribution of the odd part. The case 
of (5 —> 00 is equal to the Neumann boundary condition that separates the system in to the 
two regions, positive and negative part of the real line. Of course having two solutions is an 
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indicator of degeneracy coming from the parity invariancy of the system in this hmit. 5 = 
is the free particle case and one can see that only the even parity has contribution in the area 
distribution. 



5 Conclusion and discussion 



In this paper we found the area distribution of the elastic Brownian motion in some limits. 
Our method was based on the equality of this process with the self adjoint Hamiltonian of the 
quantum particle on the half line. The corresponding Hamiltonian for the area distribution 
is just the Hamiltonian with linear potential. The eigenvalues of this Hamiltonian after self 
adjoint extension satisfy a transcendental equation and so for the generic case the distribution 
of the area is not available, however, in some limits the calculation is tractable. We found 
perturbatively the area distribution of the Brownian excursion and the Brownian meander 
in the presence of the weekly reflecting barrier. By using self adjoint extension we found a 
unified method to classify different possible area distribution for the Brownian motion in the 
presence of the boundary. Some possible applications in diffusion phenomena in the presence 
of disorder and interacting interfaces were also discussed. 

We did similar calculations for the Brownian motion on the pointed line. The self ad- 
joint Hamiltonian in this case has three independent parameters and the eigenvalues of the 
Hamiltonian satisfy the same transcendental equation in some interesting limits. Similar 
calculations can be useful in describing different distributions of diffusing particles in the 
presence of point disorder. 

There are plenty of questions remained to be answered in the study of the elastic Brow- 
nian motion by using quantum mechanical technics. Some of them are: The case of the 
maximal height of p non-intersecting Brownian excursions and Brownian bridges is also in- 
teresting to be calculated [19], the possible connection of this study to the interacting domain 
walls of elementary topological excitations in the commensurate adsorbed phases close to the 
commensurate-incommensurate transition can be interesting. The other example is the dis- 
tribution of the time to reach the maximum [21] . Unfortunately the eigenvalue equations for 
the above cases are transcendental as we faced for one example in the end of section 3 and so 
it is impossible to get a closed formula for the distributions, however, the exact calculations 
in some limits are may be possible. The other example is the distribution of the time spent 
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by the particle on the positive side of the origin out of the total time t. This distribution 
was first calculated by Levy in the case of the Brownian motion [22j. For the pointed line 
the equations are again transcendental and need to be solved by numerical calculation. We 
mostly focused on the distributions in one dimension however one can also try to calculate 
the different distributions like the algebraic area distribution, winding number distribution 
in the pointed two dimensional space, the number of defects could be finite or infinite I^. 
We believe that the method of self adjoint extension in quantum mechanics can be useful 
to calculate such kind of distributions. It is also interesting to check our results with the 
numerical calculations. 
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